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ABSTRACT 

Aims. The aim of the present work is to study the overall mass distribution of the galaxy clusters Abell 1351 and Abell 1995 
using weak gravitational lensing. These clusters have got a very different mass structure and dynamical state, and are the two 
extremes from a larger sample of 38 X-ray luminous clusters of similar size and redshift. 

Methods. We measure shear values of faint background galaxies and correct for PSF anisotropies using the KSB+ method. 
Two-dimensional mass maps of the clusters are created using a finite-field mass reconstruction algorithm, and verified with 
aperture mass statistics. The masses inferred from the reconstructions are compared to those obtained from fitting spherically 
symmetric SIS- and NFW-models to the tangential shear profiles. We discuss the NFW concentration parameters in detail. 
Results. From the mass reconstructions we infer M2d(< r2oo)-masses of 11.7 ± 3.1 x lO 14 ftf o 1 M and 10.5 ± 2.7 x W 14 h^M Q 
for Abell 1351 and Abell 1995, respectively. About 3' north-east of the main mass peak of Abell 1351 we detect a significant 
secondary peak in the mass reconstruction as well as by aperture mass statistics. This feature is also traced by clusters members 
selected by means of their V — I colour, and hence is likely a real sub-structure of Abell 1351. From our fits to the tangential 
shear we infer masses on the order of M 200 ~ 8 - 9 x W 14 h^M Q (Abell 1351) and M 200 ~ 5 - 6 x lO 14 fef o 1 M (Abell 1995). 
The concentration parameters remain poorly constrained by our weak lensing analysis. 

Key words. Gravitational lensing - Cosmology: dark matter - Galaxies: clusters: individual: Abell 1351 and Abell 1995 



1. Introduction 

Comprising the most massive gravitationally bound struc- 
tures in the Universe, galaxy clusters are essential in 
providing a deeper understanding of the properties of 
dark matter. The recent pa pers covering IE 06 57 — 
558 (the Bulle t Cluster) by Iciowe et all l|2006al ) and 
Abell 520 by iMahdavi et all (|2007t l demonstrate the 
importance of gravitational lensing for understanding 
the matter content of our Universe. Lensing studies 
of galaxy clusters provide a powerful way to identify 
large density peaks in the Universe, independent of 
their b aryonic content dMaturi et a .112007 : Mivazaki et alJ 



Wittman et alJ l200fi IClowe et all l2006bl : iDahle et al 
20031 : ISchneidenll996l l. 



2007t iGavazzi fc Soucaill l2007t ISchirmer et all l2007t 



Send offprint requests to: K. Holhjem 

* Based on observations obtained with CFH12K at the 
Canada-France-Hawaii Telescope (CFHT) which is operated 
by the National Research Council of Canada, the Institut des 
Sciences de l'Univers of the Centre National de la Recherche 
Scientifique and the University of Hawaii. 



Neither the nature nor the dynamical state of the grav- 
itating matter affect the mass estimates obtained through 
gravitational lensing. These mass measurements are only 
changed by gravitation and the geometrical configuration 
between observer, lens, and source. Although this makes 
lensing a unique tool, such measurements can be biased by 
e.g. different mass concentrations along the line-of-sight. 

In this paper we analyse two clusters of galaxies at in- 
termediate redshift (z = 0.32), Abell 1351 and Abell 1995. 
These clusters are selected from a l ensing study of 3 8 
highly X-ray luminous galaxy clusters 1 Dahle et al. 2002h . 
They were chosen for further investigation as they rep- 
resent the two extremes in this cluster sample, regard- 
ing mass distribution and dynamical state. Re-observing 
the clusters with the wide-field camera CFH12K at 
the Canada-France-Hawaii Telescope (CFHT) provided a 
larger field of view (42f2 x 28fl) than that employed by 
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Dahle et all <|2002h (18.'8 x 18.'8). This allows us to map 
the clusters to larger radii than pr eviously possi ble 



equals 0'.'95 and l'/15 for Abell 1351 and Abell 1995. The 
number density of the lensed background galaxies is 16 
The KSB+ m ethod feaiser et alJ 19951 arcmin -2 for both clusters, and the ellipticity dispersion 



Luppino fc Kaiserl 119971 : iHoekstra et al.lll998l ) is used to (after PSF correction) for each component is a g = 0.43 

and a g = 0.51 for Abell 1351 and Abell 1995, respec- 
tively. The larger dispersion for Abell 1995 is explained 
by the 20% larger image seeing, which enlarges the PSF 
correction factors and their uncertainties. 

The CFH12K mosaic camera covers a field of 12k x 8k 
pixels in total, representing an area of 42.'2 x 28.'1 on the 
sky. The pixel scale is 0'.'206 when mounted at the CFHT 
prime focus. 

In add ition we make use of the K-band data ob- 
tained by Dahle et al. ( 2002h in order to verify neigh- 
bouring peaks present in our two-dimensional mass maps 
described in Sect. HI These data were obtained at the 
2.24 m University of Hawaii Telescope using the UH8K 
mosaic camera, covering an area of 4k x 4k pixels (re- 
binned 2x2) mapping 18'8 x 18'8 of the sky. Each im- 
age has a total exposure time of 12 600 s, resulting in an 
image depth comparable to our /-band data, with 5<r lim- 
iting magnitudes of V = 25.9 ± 0.1 and V = 25.8 ± 0.1 for 
Abell 1351 and Abell 1995, respectively. Further details 
on the reduction pro cess and co-additio n of the K-band 
data can be found in iDahle et alJ l|2002( l . 

As Abell 1351 and Abell 1995 are both located at red- 
shift z — 0.32, they have a similar correspondence between 
physical and angular scale, given as IMpc = 215". 



recover the shear values of faint background galaxies in 
the images. Using a fin ite-field reconstruction technique 
i Seitz fc Schneiderll200ll . SS01) we derive two-dimensional 
mass maps, visualising the surface mass distributions of 
Abell 1351 and Abell 1995. We also apply aperture m ass 
statistics (M ap ; Schneider 19961 : Schirmer et al. 2007 ) to 
our data, comparing the results to confirm mass peak 
detections. Finally, by fitting predicted shear values from 
theoretical models to the shapes of the lensed galaxies 
we estimate the cluster masses. We assume the clusters 
to be spherically symmetric and to have density profiles 
following either a singular isotherma l sphere (SIS) or a 
Navarro. Frenk. fc Whitel <|l997t Il995l : NFW) model. 



Data processing and analysis are carried out using 
mainly the IMCAT software packag^l, Kaiser's July 2005 
version for Macintosh. IMCAT is a tool specially designed 
for weak lensing purposes, and is optimised for shape mea- 
surements of faint galaxies. It processes both FITS files and 
object catalogues. 

The outline of this paper is as follows. In Sect. [2] we 
summarise the observations and software used for our 
study, as well as the data reduction. In Sect. [3] we de- 
scribe the shear reconstruction, focusing on shape esti- 
mates and Point Spread Function (PSF) corrections. In 
Sect. H] we present the clusters' surface mass density maps 
and attempt to verify the detected mass peaks using V — I 
colours and M ap . By comparing the measured shear pro- 
files to theoretical expectations, we model the lensing data 
in Sect. HJ Finally we present and discuss our results in 
Sect. El and our conclusions in Sect. El 

Throughout this paper we assume a ACDM 
cosmology, with Qm = 0.3, — 0.7, and 

hjo = Hq/(70 km s _1 Mpc -1 ). Errors are given on 
the la level. 

2. Observations and data reduction 

2.1. Data acquisition 

The galaxy clusters Abell 1351 and Abell 1995 are 
centred at the positions ll h 42 ro 30?7 + 58°32'21" and 
14 h 52 m 50?4 + 58°02'48", respectively. They were ob- 
served with the 3.6 m Canada-France-Hawaii Telescope 
(CFHT) on the 4 nights of 7-11 May, 2000, using the wide- 
field CCD mosaic camera CFH12K. A total exposure time 
of 5400s were obtained for both clusters, all in the /-band 
filter. However, due to seeing > 1" three exposures were 
rejected from the Abell 1351 data set, resulting in total 
exposure times of 3600 s and 5400 s for Abell 1351 and 
Abell 1995, respectively. This corresponds to a 5er limiting 
magnitude of / = 25.2±0.1 for point sources in both point- 
ings. The seeing in the final co-added images respectively 



1 IMCAT is developed by Nick Kaiser (kaiser@hawaii.edu), 
http : //www. if a.hawaii . edu/~kaiser/ imcat/ 



2.2. Image processing 
2.2.1. Pre-processing 

To remove the bias level in each frame we used the mean 
value of the overscan region from the corresponding chip. 
The flat-fielding was carried out using a master night time 
flat, made from averaging 56 night time exposures; most 
of them so-called "blank" fields and all well displaced from 
each other. The fringing that occurs in /-band exposures 
is also captured in this type of flat, and upon dividing the 
object exposures by it the fringes were cleanly removed. To 
estimate the background level in the exposures, we used 
the heights of the minima of the sky level present to create 
a model for each individual frame. After subtraction the 
median sky level was set to zero. 

As fringing is an additive effect and not a multiplica- 
tive one, ideally the fringes should be subtracted. Since we 
had no twilight flats available, standard defringing could 
not be performed. The photometric error introduced by 
division is negligible, as the amplitude of the fringes com- 
pared to the sky background after flat-fielding was on the 
order of 2%. However, since fringing acts mostly on small 
angular scales, its treatment will affect the shapes of the 
small and faint background galaxies used for weak lensing. 
To investigate this we obtained a set of 10 ar chival images 
of the Deep3 field I Hildebrandt et al. 2006), taken with 
the Wide Field Imager at the 2.2m MPG/ESO-telescope 
through their /-band filter. As the Deep3 field does not 
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contain any massive clusters it is very well suited for this 
test. Two different coadded images were created. In the 
first case the data were flat-fielded using twilight flats, 
and then a fringing model was created from the flat- 
fielded data and subtracted. The second coadded image 
was processed in the same way as our CFHT data, i.e. 
the data were flat-fielded and fringe-corrected by divi- 
sion of a night-time flat. We then measured the shapes 
of a common set of ~ 12000 galaxies in both images (see 
Sect. [3]) and created two mass reconstructions, using the 
same technique and smoothing scale as for Abell 1351 and 
Abell 1995 (see Sect.|3|. We find that the rms of the differ- 
ence of the two mass maps is a factor of 2.5 smaller than 
the noise of the individual mass maps, mainly caused by 
the intrinsic ellipticities of background galaxies. The effect 
in our CFH12K data is much smaller, as the CFH12K I- 
band filter has a blue cut-on at around 730 nm and a 
red cut-off at 950nm. The ESO /-band filter on the other 
hand opens at 800nm and has no cut-off on the red side. 
Hence the fringing amplitude in the comparison data set 
from ESO is up to 5 times higher than in CFHT. We con- 
clude that our analysis of Abell 1351 and Abell 1995 is 
not affected by our fringe correction. 

2.2.2. Masking 

The CFH12K mosaic contains some bad pixels and 
columns, in particular two of the CCDs suffer from this. 
By using Nick Kaiser's ready made CFH12K mask^l as 
global masks, all bad areas were ensured to be ignored. 
An additional patch of 219 bad columns in CCD00 was 
also added to the global masks. We did not make further 
individual masks for each exposure, as most spurious de- 
tections are filtered out during the astrometric calibration. 
In addition, suspicious objects are rejected from the final 
object catalogue by visual examination. 

2.2.3. Astrometric calibration 

Wide-field data typically do not have a simple relation 
between the sky coordinates and those of the detector. A 
mapping from pixel coordinates onto a planar projection 
of the sky needs therefore be performed. We solved for this 
through a series of steps. 

First, all objects in each exposure were detected and 
aperture photometry carried out. By plotting r g (~half- 
light radius) vs. instrumental magnitude of the objects 
in each exposure, we extract the moderately bright, non- 
saturated stars suitable for deriving an astrometric so- 
lution. Their weig hted ellipticity parameters (defined by 



Kaiser et al. 19951 ). ei : 2, are in addition plotted and eye- 



balled, selecting the main clustering of objects to ensure 
a catalogue containing purely stars (see Fig. Q] for an ex- 
ample) . 

Left with star catalogues we now compute the trans- 
formation parameters needed using information from the 



USNO-B1.0 catalogue (|Monet et al.ll2003l l. However, as 
many of the USNO-B1.0 stars are saturated in our im- 
ages, we extended our reference catalogue by detecting 
more sources in FITS files derived from the Digitized Sky 
Survey (DSSfl We match the target catalogues to the 
reference catalogue, then solve for a set of low-order spa- 
tial polynomials mapping the images onto each other, by 
repeatedly refining the least squares minimisations using 
outlier rejection. 



2.2.4. Final master image and object catalogue 

Gain and quantum efficiency variations between chips and 
differential extinction between exposures are solved for 
by least squares minimisations. The extinction corrections 
between the exposures were very small, typically ~ 0.01 
mag, whereas the zero-point offsets between the chips were 
~ 0.1 mag. As an accurate absolute photometric calibra- 
tion is not necessary for the present work, we adopted 
standard Landolt magnitude zeropoints. 

The co-addition is done after magnitude corrections 
are applied to the data. In addition cosmic rays are masked 
out, before the median image is computed and the back- 
ground flattened. A master object catalogue is created, 
where each object's WCS coordinates are calculated from 
the astrometric solution, and their ellipticity parameters 
computed (see Sect. [3]). Finally we mask out false detec- 
tions by over-plotting the objects onto the image, hence 
ensuring a final object catalogue free from spurious detec- 
tions. 



3. Shear measurements 

Identifying weak lensing effects requires measuring the el- 
lipticities of a large number of faint background galaxies. 
The main source of noise in weak lensing analysis is the 
intrinsic ellipticities of these galaxies. To distinguish be- 
tween distorted images resulting from a weak lens and the 
usual distribution of shapes existing in an unlensed galaxy 
population, the ellipticities are examined for a systematic 
change. In particular, a tangential alignment of the galaxy 
shapes around the cluster centre would confirm the exis- 
tence of a weak lens. 

An additional source of error comes from the faint 
background galaxies being smeared by the PSF, caused 
by atmospheric turbulence and optical aberrations. The 
weak shear signal is hence diluted because this smearing 
will cause the galaxy images to appear more circular than 
before the smearing. In addition, PSF anisotropies distort 
the images, causing the galaxies to appear more ellipti- 
cal, hence introducing false shear signals. It is crucial that 
these PSF effects are corrected for. 



http : //www. if a.hawaii . edu/~kaiser/cfhl2k/masks/ 3 http : //archive . stsci . edu/dss/ index.html 
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Fig. 1: Example of size-magnitude (left) and weighted ellipticity parameter diagram (right) for an arbitrary exposure. The 
moderately bright, non-saturated stars are chosen from within the squares and utilised in the astrometric solution, the right 
plot containing only the stars chosen from the left plot. As square is the only option IMCAT offers in selecting objects, this is 
also applied to the ellipticity diagram. Although a circle might yield better results, the difference is considered negligible. The 
stars utilised in the PSF corrections (Sect. [3~Tjl are selected the same way. 



3.1. PSF corrections 



Follow in g the KSB+ meth o d dev eloped by iKaiser et al 



Luppino &; Kaiserl l|1997f ). and iHoekstra et al 



1998), we present a short summary of our implementa- 
tion below. KSB+ inverts the effects of PSF smearing and 
anisotropy on objects in an image, presenting a method 
to recover the true shear. 

Ignoring the effects of photon noise, it is possible to 
express the observed ellipticity of a galaxy as 



Pi 



(1) 



where the first term represents the intrinsic ellipticity of 
the galaxy, the second term the shift in ellipticity caused 
by gravitational shear, and the third term the smear- 
ing of the galaxy imag e from the anisotropic PSF (see 
Luppino fe Kaiserl 1997 with additional corrections from 
Hoe kstra et al J11998J for a thorough deduction of this equa- 
tion). Present in eq. (HJ are the pre-seeing shear polaris- 
ability tensor, -PTg, and the smear polarisability tensor, 
P^a- The latter is calculated for each object together with 
e° bs , the post-seeing shear polarisability tensor P^g, and 
the centroid of the object. 

Because stars are foreground objects (gp = 0) and in- 
trinsically circular (e s a = 0), applying eq. Jl} to stellar 
objects provides a measure of the total PSF anisotropy, 
Pp. This is calculated from bright stars selected from the 
final object catalogue (Sect. f2T2T4|) . The PSF corrections 
are then calculated for all individual objects and correc- 
tions applied respectively. 

The ellipticities of the stars were fitted to a sixth- 
order Taylor series expansion. When comparing mass and 
B-mode maps (Sect. |4]) for fits of different orders, there 



was little change with the order of fit. Over the whole 
field, 410 and 530 stars were used in the fitting process for 
Abell 1351 and Abell 1995 respectively. Fig. [2] shows the 
ellipticities of the stars before and after PSF corrections. 

The pre-seeing shear polarisability tensor, Pj^, is de- 
fined in KSB+ to be 



p7 _ psh 
r aP — r afi 



psm 



(P S 



t \ X psh* 



(2) 



where the asterisk denotes P^% and P™ 



^5 1 50 

applied to stel- 
lar objects. From eq. Jl]) we see that the reduced shear, 
.9 = 7/( 1 -«0) is g iven b Y 

obs 



/ n7 \ — 1 f obs nsm„ 1 



(3) 

Following the approach by Wold et al. 1 20021. we as- 
sume the PSF is close to circular after the correc- 
tion, and the polarisabilities can be approximated by 
P af} = |(P U + P 22 )5 aP . The average of p sh */p sm * can be 
calculated as 



P 



Sli+ 



1 



iVs, 



E 



psb 



sh* 



psh* 
r 22 



P 



11 



psm* 
r 22 



(4) 



where we use the median value rather than the mean to 
minimise the effect of outliers. Hence eq. |(2]) turns into 

i / psh* 



p-f = + p*) _ l ( p i7 + p-) ^ 



(5) 



As P^p and P^™ are already calculated, this equation 
is easily solved, and eq. |(3]) gives us an estimate of the 
gravitational shear of each object. [Hoekstra et al. (1998) 
show that estimating the PSF dilution for each individ- 
ual galaxy introduces additional noise. We therefore fol- 
low their approach by determining P 7 as a function of 
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Fig. 2: Ellipticities of the stars in the field of Abell 1351 (top) and Abell 1995 (bottom) before and after corrections for PSF 
anisotropies. The stars initially have systematic ellipticities up to ~ 7 — 9% in one direction. The PSF corrections reduce these 
effects to typically < 1.5%. 



magnitude and galaxy size. We determine the median P 7 
within 15 bins in an r g -magnitude diagram, where the cen- 
tral 4 bins contain ~ 4000 galaxies /bin and the outer ones 
~ 200 galaxies/bin. We then compute one correction fac- 
tor for each bin using eq. |(5]) , and apply this to all galaxies 
within the corresponding bin. 



The faintest and smallest galaxies are more affected 
by seeing than the larger galaxies, giving them a poorer 
shape determination and a larger correction factor. Such 
galaxies are therefore of less importance. To account for 
this, a normalised weight, 
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is calculated for each bin i and assigned to the correspond- 
ing galaxies. Here, a ei is the observed dispersion of galaxy 
ellipticities. 

4. Mass reconstruction 

We select background galaxies with 6 < S/N < 100 for 
the creation of our mass maps. These are reconstructed 
from the galaxies' s hapes using the fin i te-fiel d inversion 
method presented bv lSeitz fc Schneider 1 200lL SS01). For 
this method a smoothed shear field is calculated on a grid 
using a modified Gaussian filter. The algorithm then it- 
eratively computes a quantity K{9) := ln[l — k(9)] which 
is determined up to an additive constant due to the mass 
sheet degeneracy. We break the latter by assuming that 
the average convergence vanishes along the border of the 
wide field of view. The width of the Gaussian term in 
the filter is set to 2.'0, resulting in an effective smoothing 
length of about 1'.6. 

In order to evaluate the noise of the mass maps, we 
computed 2000 mass maps for each cluster based on ran- 
domised galaxy orientations, keeping their positions fixed. 
As the cluster lens signal increases the ellipticities of galax- 
ies, this would lead to an overestimation of the noise at the 
cluster position. We roughly correct for this effect by sub- 
tracting the expected SIS tangential shear signal, deter- 
mined from the clusters' known velocity dispersions (see 
Table[l]). Since the singularity of the SIS can lead to overly 
large corrections close to the cluster centre, we limited the 
maximum correction factor allowed to 0.5 in each elliptic- 
ity component. This affected less than 5 galaxies in both 
fields. The true mass maps were then divided by the noise 
maps obtained from the randomised mass maps to create 
the S/N-maps seen in Fig. [31 

Abell 1351 and Abell 1995 are detected with 
a S/N of 5.3 and 5.2, respectively. Upon integrat- 



ing the k maps within r 2 no 



1.69 h 



TO 



Mpc 



(1.50 /ifo 1 Mpc) for Abell 1351 (Abell 1995), we find 
total masses of M 2D (< r 200 ) = 11.7 ± 3.1 x 10 14 h^M & 
and M 2D (< r 2Q0 ) = 10.5 ± 2.7 x 10 14 h^M Q for the clus- 
ters, respectively. The r 2 no radii are taken from what we 
consider to be the best NFW fits to the data (see Table Q] 
and Sect. 15.21 ). while the errors are determined from inte- 
grating the same areas in the 2000 noise maps. 

In Fig. [3] the B-modes in both cluster fields are also 
shown, computed by re peating the k reconstru ction with 
all galaxies rotated 45° I Crittenden et al. 20021 ). Provided 
that the lensing data are free from systematics and that 
the noise (intrinsic ellipticities) is Gaussian, this B-mode 
map should be consistent with Gaussian noise. Given the 
effective filter scale of 1.'6, we can place about 380 indepen- 
dent peaks in the CFH12K field. Thus one would expect 
1.1 noise peaks above 3<r in the field. A more realistic es- 
timate comes from the 2000 randomisations as these are 
based on the real ellipticity and spatial distribution. We 
expect 1.4 (1.6) such peaks for Abell 1351 (Abell 1995). 
In the real B-mode maps we find 3 peaks for each of the 
clusters. This is insignificant, as in our randomisations at 



least 3 such peaks appear per field in 20% of the cases. 
In case of Abell 1995 the largest B-mode peak has a sig- 
nificance of 3.9ct. Its B-modes appear generally somewhat 
larger than for Abell 1351, which has no B-mode peaks 
higher than 3.5o\ 

4.1. Mass and galaxy density distributions 

In order to compare the surface density maps with the dis- 
tribution of_£husJej_gahMie^we extract the red sequence 
(see e.g. iGladders fe Yeel [2000) and investigate the dis- 
tribution of the galaxies selected. We combine the /-b and 
data with the V-band images from lDahle et al. 120021). To 
match the two data sets we re-sampled both images to the 
same pixel scale, resulting in a common area of 18f 5 x 18.'5 
on the sky. The F-band image seeing was around Of! 7, and 
thus consistently better than in the /-band. The F-band 
data were therefore convolved to match the seeing in the 
/-band data, 0'.'95 for Abell 1351 and l'/15 for Abell 1995. 
Aperture photometry w as carried out using SExtractor 
( Bertin fe Arnoutsll 19961 ) in double- image mode. The deep 
/-band images served as the detection image, providing us 
with a target list with defined coordinates. At these posi- 
tions we integrated the flux in a 3" wide aperture in each 
of the V- and /-band images. Plotting the galaxies in a 
colour-magnitude diagram will then in principle provide 
enough information to separate the red early-type cluster 
members from the other galaxies. 

Each cluster's red sequence does not clearly stand out 
from the V — I vs. / diagram when all objects are plot- 
ted. We therefore select only galaxies within a radius of 
3' of the brightest cluster galaxy (BCG) for the colour- 
magnitude diagram (see Fig. HJ, and detect the red se- 
quence at 1.4 < V — I < 1.9 for both clusters. The selec- 
tion criteria indicated by the box in each plot are then 
applied to the entire object catalogue. The number den- 
sity of the galaxies selected is then calculated as a function 
of position and over-plotted onto the central 9 arcmin of 
the mass maps (see Fig. [5]). 

We normalise the number density maps by the fluc- 
tuation measured in the field outside the clusters. The 
centres of Abell 1351 and Abell 1995 are then detected 
with 15.5cr and 9.9<r significance, respectively. The posi- 
tional offsets between mass centres, BCGs, and centres of 
galaxy density distributions are in the range of 30" — 55" 
for both clusters, and are due to noise in the mass maps. 
Changing the width of the Gaussian kernel in the finite- 
field reconstruction algorithm shows that the peak centres 
can drift by up to 20" from the mean position. These off- 
se ts are consist e nt with other results in the literature, such 



as iClowe et alJ (|2006al ) , who observed offsets on the order 
of 10" between the lensing peaks (of higher S/N than ours) 
and the optical centres of the bullet cluster. Positional off- 
sets of 50" are common in t he sa mple of 70 shear-selected 
clusters bv lSchirmer et al.1 |2007j). 

The cluster galaxy distribution resembles well the mass 
distribution in the central part of Abell 1351. It extends 
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Fig. 3: The projected surface mass densities and B-modes for both clusters in the full CFH12K field, using a finite-field mass 
reconstruction. The maps show the S/N of the clusters, with contours starting at Oct and increasing in steps of la. The clusters 
are detected at significance levels of 5.3ct (Abell 1351) and 5.2a (Abell 1995). The smaller peaks A, B, and C have S/N-ratios 
of 4.2, 3.8, and 3.8, respectively. The effective smoothing length of the reconstructions is 1'.6. The highest B-mode peak of 
Abell 1351 (Abell 1995) has a S/N-ratio of 3.5 (3.9). 



significantly towards the north-east, a feature also seen 
in the mass map where we find a local maximum which 
we refer to as peak A (see Sect. 14. 2p . The galaxy distri- 
bution of Abell 1995 appears elliptical and elongated in 
the northeast-southwest direction. This property is not 
reflected in the mass map where the peak is of rather cir- 
cular appearance. 

To check the integrity of our mass reconstructions fur- 
the r, we applied t h e pea k finder (^-statistics) developed 
by ISchirmer et al 1 2007 ). This method detects areas of 
enha nced tangential shear using the aperture mass statis- 
tics l|Schneiderlll996t ). Since it uses a filter function that 
mimics the tangential shear profile of galaxy clusters it is 
well suited for detecting mass concentrations. With this 
method we recover Abell 1351 at the 7.0a level in a 10' 
wide filter, and Abell 1 995 with 6.1a for a 7 ' filter. The 
filter shape parameter ijSchirmer et al.l 120071 ) was chosen 
as x c = 0.2 in both cases. We find that the S-statistics is 
elongated in the same way as the mass reconstruction for 



Abell 1351, extending towards peak A. We evaluate the 
significance of this possible sub-structure in the following. 

4.2. Smaller mass peaks in the fields 

In the mass reconstructions two neighbouring peaks A and 
B are detected around Abell 1351, and another one (peak 
C) in the field of Abell 1995. Their S/N-ratios are 4.2, 
3.8, and 3.8, respectively. Using the 2000 noise randomi- 
sations for each field we find that the probability of a noise 
peak larger than 4.2a (3.8a) in the field of Abell 1351 is 
0.8% (6.8%), respectively. The corresponding probability 
for peak C in the Abell 1995 data is 7.0%. These are some- 
what larger than what would be expected from purely 
Gaussian noise. 

Hence the only significant sub-structure we detect 
in the mass reconstruction is peak A near Abell 1351. 
Looking at the contours in the upper right panel of Fig. [5] 
one can see that the S'-statistics trace this structure as 
well at the 4.5 — 5.0a-level. We note that we recover this 
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Fig. 4: Colour-magnitude diagrams for Abell 1351 (left) and Abell 1995 (right), where only galaxies within 3' from the BCG are 
plotted. The red sequence can be seen for 1.4 < V — I < 1.9 for both clusters, the box indicating our selection criteria. 



sub-structure over a broad range of filter scales in the S- 
statistics and hence think that it is a real feature in the 
mass distribution of Abell 1351. 

Out of the broad range (1' — 15') of filters probed with 
the .S'-statistics, peak B is detected only once with 4. Oct in 
the 4' wide filter for x c = 0.5. It h as the typical ch a racter - 
istics of the dark peaks found by ISchirmer et alJ (|2007n , 
i.e. it is not associated with any over-density of galaxies. 
Hence it is most likely a noise peak. 

In the Abell 1995 field we could not detect any other 
peaks using the S'-statistics. Since the B-modes for those 
data show a maximum of 3.9ct near peak C (at 3.8ct), we 
consider it to be a noise peak. As it also lies outside the 
area covered by the F-band, we could not check for over- 
densities of red galaxies at this position. 

5. Modelling the lensing data 

Comparing observed distortions in the background galax- 
ies to those predicted by theoretical density profiles en- 
ables us to estimate the mass of a galaxy cluster. Using 
^-minimisations of SIS and NFW models, we first deter- 
mine the best fit parameters and then calculate the cluster 
masses. 

The theoretical profiles are both spherically symmet- 
ric. We therefore average the tangential reduced shear, 
9t = 7t/(l — «0 (for r > 6>e, where 6*e is the Einstein ra- 
dius), in 17 radial bins around the cluster centre. The 
bins are logarithmically spaced, covering the entire field of 
view, and starting at r min = 150" to avoid the large con- 
tamination from cluster galaxies close to the centre of the 
field (see also Sect. I5.1| ). To determine the cluster centre, 



we tested three different positions. First we adopt the peak 
location in the mass reconstructions generated (Sect. H]). 
These coincide with the centres of the S'-statistics. Second, 
the BCG serves as cluster centre, and third we try the cen- 
tre of the galaxy density of each cluster. As the latter coin- 
cide with the BCG for Abell 1995, only two positions were 
tested for this cluster. We also considered strong lensing 
features, but found that they do not offer further insight in 
this respect (see Sect. l6.3TTl for details). We calculate (gt)i 
for each radial bin i and compare them to the theoretical 
values at the average radius of each bin, (r) -. 

When calculating the mass of a cluster, the relative dis- 
tance of the background galaxies and the lensing cluster 
is required. As we have no specific information about the 
redshifts of the background galaxies, the distances have to 
be estimated statistically. By using the photometric red- 
shift distribution of corresponding f aint galaxies from the 



Hubb le Deep Field (HDF) North l|Fernandez-Soto et al 
Il999n . we can estimate the average (3 = Da s /D s , where 
.Dds is the angular diameter distance between the lens and 
the source and D s between the observer and the source. 
The empirical relation 



</3> 



-1.21z d + 0.91 



(7) 



is deriv ed for a ACDM cosm ology analogously to the equa- 
tion of Dahle et al. ( 2002T) for an Einstein-de Sitter uni- 
verse, which also accounts for the redshift distribution of 
the source galaxies. Here denotes the redshift of the 
lens. 
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Fig. 5: The black contours show the mass reconstruction (left) and S-statistics (right) for Abell 1351 (top) and Abell 1995 
(bottom). The contours start at the 3. Oct level and increase in steps of 0.5ct. Over-plotted in white contours (red in online 
version) are the number densities of galaxies selected with the red sequence method, normalised by the rms fluctuation in the 
field outside the clusters. These contours are isodensity contours, starting at 3ct going in steps of 3ct, and peak at 15 and 10a, 
respectively. 



5.1. Cluster contamination and magnification depletion 

At small projected radii from the cluster centre our faint 
galaxy sample will contain cluster galaxies in addition to 
background galaxies. We could not discriminate faint clus- 
ter members from lensed field galaxies using V — I colours, 



hence the sample of presumed lensed background galaxies 
remains contaminated. This leads to a systematic bias of 
the shear measurements towards smaller values. 

In order to quantify this contamination, we determine 
the over-density of galaxies in the background catalogue 
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at the cluster position compared to the mean density in 
the field (an example is shown in Fig.[6]for Abell 1351). A 
contamination rate of 50% is found for both cluster cen- 
tres, vanishing for radii larger than about 4' — 5'. It is in 
fact even higher, as magnification depletion leads to a re- 
duced number density of lensed galaxies in the /-band near 
the cluster centre. However, this effect can be neglected in 
our case. From the smoothed convergence (see Sect.H]) and 
reduced shear fields we estimate the magnification using 
/i = ((1 — k)(1 — g))~ 2 . We find very similar magnifications 
for both clusters, being 1.65 at the centre and becoming in- 
distinguishable from the noise (ct m ~ 0.15) for radii larger 
than ~ 3'. The depletion of the number density of lensed 
galaxies is oc jn 2 ' 55-1 , with s = 0.15 in red filters (see e.g. 
Naravan 8z Bartelmann 199ot . At the cluster centres the 
number densities are thus reduced by a factor of ~ 0.73, 
and at a radius of 1.'5 magnification depletion becomes 
indistinguishable from the natural fluctuations in the dis- 
tribution of field galaxies. Magnification depletion hence 
only affects the innermost ~ 0.3 Mpc (100") of the clus- 
ters and can be neglected since we compare the tangential 
shear profiles to models only for radii larger than 0.5 Mpc 
(see Fig. [7|). 

In order to correct for the contamination by cluster 
galaxies, we modify the theoretical shear values. The rea- 
son for adjusting the theoretical values rather than the 
measured values is that this method is considerably eas- 
ier to implement. The correction factors are determined 
in radial bins of logarithmic spacing. One correction value 
is then calculated for each of the 17 bins in which (gt) is 
measured. By assuming the edges of the field to be ap- 
proximately free from cluster galaxies, the outermost cor- 
rection factor can be set to 1 to mimic contamination-free 
boundaries of the field. Finally the best fit is found using 
^-minimisations. 



5.2. Fitting the SIS and the NFW profiles 

Once the cluster centre is determined, the only free param- 
eter of the SIS profile is the velocity dispersion, cr„. The 
best fit of the SIS profile is determined by x 2 -minimisation 
for a range of o v values, the results being shown in Table[TJ 
The mass estimate, Mgis, for this profile is calculated at 
7*200 (the radius inside which the mean mass density of the 
cluster equals 200p cr u) found in the NFW fitting with two 
free parameters utilising the same cluster centre. 

The NFW profile is derived from fitting the density 
profiles of numerically simulated cold dark matter halos. 
It appears to give a very good description of the radial 
mass distribution inside the virial radius of a galaxy clus- 
ter. For a thorough introduction to the gravitational lens- 
ing propertie s of the NFW mass densit y profile we refer 
the reader to Wright &: Brainerdl 1 2000l ). The theoretical 
7t and K ca n be calculated analytically for the NFW den- 
sity profile 1 Bartelmann fc SchneideifeOQl ). We derive the 
best fit parameters for different values of the concentration 
parameter, c, ranging from 0.1 to 24.9 in steps of 0.1. 



Cluster galaxy contamination in faint galaxy sample 




5 10 

radius [arcmin] 

Fig. 6: Percentage of cluster galaxies in the faint galaxy cata- 
logue of Abell 1351 (that of Abell 1995 is very similar). Because 
the projected density of cluster galaxies is assumed to equal 
zero at the edge of the field, the cluster galaxy contamination 
is set to zero here by subtracting the median value outside the 
central area of the image (the large field of view makes this a 
well- working approximation). Our fitting procedure starts at 
r m i n = 150" to avoid the large cluster galaxy contamination at 
the centre. 

With the cluster centre fixed, the NFW profile has two 
free parameters, r2oo and c. We fitted our shear measure- 
ments to this profile twice; first by keeping c fixed and 
varying only r2oo to find our best fit, and second by vary- 
ing both parameters. The best fit parameters were deter- 
mined by minimising \ 2 m both cases. Based on N-b ody 
simulations of dark matter halos . iBullock et al.l (|200lh de- 
rive relations for the mean value of c as a function of red- 
shift and mass for different cosmologies. For a halo of mass 



x 1O 14 M , the relation yields 
5.8 



1.194(1 + zd) 



(8) 



(where r2oo = 1.194r V i r for a ACDM cosmological model). 
As this mass is close to the mass estimates of Abell 1351 
and Abell 1995 (Sect.H]), the weak mass dependence of c 
can be disregarded. In the second fitting process both r2oo 
and c were altered, creating a grid of c, r2oo-values. The 
best fit r2oo was first determined for each value of c, then 
the best fit c was found. The results are given in Table [TJ 
We could not find an upper limit for the concentration 
parameter of Abell 1351, independent of the cluster centre 
chosen. The same holds for Abell 1995 in case the BCG 
is chosen as the centre. We discuss this in more detail in 
Sect.EZH 

As an example we display the reduced tangential shear 
as a function of radius using the SS01 n maps as clus- 
ter centre, see the left diagrams of Fig. The measured 
values of (g t ) are given together with the best fit mod- 
els of the SIS and NFW profiles. Judging from the dia- 
grams alone, the NFW profile letting both c and r2oo vary 
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seems to represent the best fit to the clusters. However, 
the x 2 /Dof values given for each model in Table Q] show 
that the differences between the models are not statisti- 
cally significant. The differences emerging from different 
cluster centres seem to have a larger impact. The right 
diagrams of Fig. [7] show the B-modes of both clusters, i.e. 
the cross-component of (gt) as a function of radius. Both 
measurements are consistent with zero. 

6. Discussion 

X-ray studies show Abell 1351 to be a galaxy cluster ex- 
hibiting significant dynamical activity and undergoing a 



major merger event l| Allen et al.ll2003l ). which indicates a 



cluster still in its formation phase. Analyses assuming a 
relaxed cluster will hence naturally differ from weak lens- 
ing analyses, where no assumption is made about the dy- 
namical st ate of the cluster. One example is the virial 
analysis by llrgens" et alJ ()2002h . where an unusually high 



velocity dispersion of a v = 1680 _ 22g km s is obtained 
for Abell 1351, based on radial velocity measurements of 
17 cluster galaxies. Such a high velocity dispersion is not 
uncommon in merging systems. If for example two smaller 
clusters with low velocity dispersions fall towards each 
other along the line of sight with a velocity comparable to 
or larger than their a v , then a very large total a v would 
be inferred, with a correspondingly overestimated virial 
mass. The cluster CL0056. 03-37.55 is a good example for 
such a system (see Schirmer et al.ll2003l ). 

Abell 1995 is, unlike Abell 1351, classified as a re- 
laxed cluster in dynamical equilibrium ( Pedersen 8z Dahlel 
2006). X-ray studies and virial analyses of this clus- 
ter are hence also more compatible with lensing studies 
i Patel et alj|2000l llrgens et alj|2002l ). The projected two- 
dimensional distribution of cluster galaxies in Abell 1995 
is clearly elliptical (see Fig. [5j) , whereas the central lensing 
mass distribution is circular. 

6.1. The mass estimates 

The mass distributions of Abell 1351 and Abell 1995 are 
estimated assuming that the clusters follow spherically 
symmetric SIS or NFW profiles. Although an elliptical 
mass pr ofile might yield more accurate cluster mass esti 
mations. 



The concentration parameters do not change significantly 
by this boosting of ellipticities. Since we do not know by 
exactly how much our shear values are underestimated 
this was merely an attempt to quantify this effect on our 
data, and is not taken into account in the results presented 
in this paper. 

Dahle et al. 1 20021 ) obtain weak lensing estimates of 



the cluster velocity dispersion of several clusters using 
an SIS model and assuming an Einstein-de Sitter 
Universe. Their results are given for Abell 1351 
as o v = 1410+^° kms" 1 and for Abell 1995 as 
a v = 1240 ± 80 km s~ , and do not agree with our 
results. However, there are several imp o rtant differences 

(|2002h and our 



in methodology between iDahle et al. 
work. As mentioned above, the assumed cosmological 
model is different. Also, these authors approximate 
9t = 7ti whereas we use <?t =7t/(l — K ) in our fits and 
mass reconstructions. Finally, the shear estimator of 



Kaiser (2000) u sed by IDahle et all (20021 ) is shown by 



Heymans et al 



response 



1 2006n to ha ve a non-lin ear 

to shear. A re-analysis of the Dahle et al. 1 20021 ) data, 
taking all these effects into account, yields new values of 
a v = 1410 ± 90 km s _1 and a v = 1000 ± 100 km s" 1 for 
Abell 1351 and Abell 1995, respectively. Hence there still 
remains a systematic discrepancy between the results for 
Abell 1351, while the measurements for Abell 1995 agree 
within error bars. 

A rema i ning difference between our 
Dahle et all (2002) is the maximum radius, 



work and 
r ma x, to which 

the shear is measured, given by the field of view of the 
detector. Changing r max in our Abell 1351 data to 550" 
(as this is the value used by Dahle et al. 20021 ) yields 



a„ = 1240 ± 105 k m s , consistent with the re-analysed 



Dahle et al. 1 20021 ) values within error bars. 



Allen et all (2003) use the IDahle et all (20021 ) obser- 



vations to obtain a weak lensing mass estimate apply- 
ing the NFW model to a ACDM cosmology. Their re- 

}Mn for Abell 1351 



suits give M 200 = 30.2+^ x 10 



14 



hj iviq 



and M 2 oo = 
values are 



lAllen et al 



14.4+^ x 10 14 h^M e for Abell 1995. These 
l igh c ompared to the results of this study. 
(|2003h use a fixed concentration parameter 
fitting process, c = 5. Applying this value to 

The discrep 



in the nttmg process, c 
our data yields minimal ch anges in M-m 



Dietrich et alJ jSj) experience that the results ancies he . n . c ^ °^ inat< : from l Allen et al. | jjjjfl3|) utilising 



from fitting a singular isothermal ellipse (SIE) model de- 
pend strongly on the initial values chosen for the minimi- 
sation routines. We therefore decided not to fit the SIE 
pr ofile to our clust ers . 

Heymans et alJ 1 20061 ) demonstrate in the Shear 
TEsting Program (STEP) that the KSB+ PSF correc- 
tion tends to systematically underestimate the shear val- 
ues ~ 10 — 15%. To measure how much this affects our 
data, we calculate an upper limit for our mass estimates 
by increasing the ellipticities with 15% and repeating the 
fitting process. We find that the underestimation of shear 
leads to an underestimation of the total cluster mass with 
a maximum of 20%, which is within the initial error bars. 



r 2 oo = 2.69 



„ l9 h^Mpc and r 200 = 2.07^ h^Mpc 
(priv. comm.) for Abell 1351 and Abell 1995, respectively, 
as these values are larger than the r2oo values we obtain 
as our best fit. 



6.2. The concentration parameter 

The mass density of a cluster with a small concentra- 
tion parameter decreases slower when going to larger radii 
than for a cluster with a large c value ( Wright fc Brainerdl 
2000). Although unconstrained upwards, we find a lower 
limit of c > 4.5 for Abell 1351. As is also seen from the 
radial dependence of the shear in Fig. [7] (top left), the 
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Fig. 7: Left: Reduced tangential shear as a function of radius for Abell 1351 (top) and Abell 1995 (bottom) using the SS01 k 
map as cluster centre (the other centres show very similar figures). The averaged gravitational lensing distortions of background 
galaxies are shown as points with ler error bars. The lines indicate the best fit models; the dashed line representing the SIS 
profile, the dot-dashed line the NFW profile keeping c fixed, and the solid line the NFW profile with two free parameters. 
It should be noted that these lines represent the theoretical values after modifications from cluster galaxy contamination are 
applied to each bin independently of each other (see Sect. [5}. Right: Cross-component of the reduced tangential shear as a 
function of radius for Abell 1351 (top) and Abell 1995 (bottom). This signal should disappear if {g t ) is caused by lensing only, 
and it is seen that the measurements are close to zero for both clusters. 



mass distribution of Abell 1351 concentrates around the 
cluster centre, indicating that its concentration parameter 
is significantly higher than that of Abell 1995. The found 
values for c of Abell 1995 (see Table [J) suggest that its 
mass is spread more evenly to larger radii, which is also 
seen in Fig. [7] (bottom left). 



Fr om their aperture mass calculations iDahle et al 
1 20021 ) find that most of the mass of Abell 1995 is 
contained within r ~ 0.9 h^Mpc (~ 200"). The mass 
of Abell 1351 shows the opposite behaviour, increasing 
evenly with radius, even at large radii. These results are 
contrary to our conclusions. As measurements at large 
radii are certain to include additional information not 



recognised close to the cluster centre, these discrepan- 
cies are likely explained by the difference in field-size 
between the two studies. By mapping the mass distri- 
bution towards a r a dius more than twice the size as 
that of Dahle et al. 1 2002h . our results are better con- 
str ained. Furth e r bias also arises from the measurements 
of IDahle et all l|20Q2t) starting from an inner radius of 
fmin = 50" (~ 0.37 h? Q Mpc), where we consider the clus- 
ter galaxy contamination to be very high, in addition to 
not including any correction for this contamination. 

Our shear values are measured from a radial cut-off, 
^min = 150", to avoid the large cluster galaxy contamina- 
tion present at small radii. Because c is estimated from 
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Table 1: Results from fitting theoretical density profiles to the measured shear values, where each sub- headline indicates the 
cluster centre around which the fitting was done. The mass estimates for the SIS profile (Msis) are calculated at the best fit 
T2oo from the NFW profile with two free parameters using the same cluster centre. The NF W profile denoted "f ixed" c refers to 
the fitting process keeping c fixed while varying only r 2 oo, where c is calculated from eq. (HJ ijBullock et al.l[200lf l. In the bottom 
we add the results from integrating the k maps within our best fit r2oo values. Note that the KSB+ PSF correction tends to 
underestimate the shear by 10 — 15%, which in turn reduces the cluster masses up to 20%. 
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the scale radius, r s = r2oo/c, it is desirable to include r s 
in the measurements (r mm < r s ) in order to obtain an ac- 
curate estimate of the concentration parameter. If this is 
not the case, c is basically unconstrained. 

This appears to be the case for Abell 1351, explain- 
ing why we were not able to derive an upper limit for its 
concentration parameter. Letting r m - ln = 150", we ensure 
a cluster galaxy contamination < 25% at this inner ra- 
dius. However, as the c parameter appears unconstrained 
under this condition, we reduce r m i n attempting to obtain 
clearer results. The problem now arising is the increasing 
contamination of cluster galaxies. Looking at Fig. [6] we 
see that at r — 120" the cluster contamination is ~ 32%, 
and at r = 100" it equals ~ 40%. Though this contamina- 
tion is accounted for during the fitting process, the con- 
tamination correction is still vulnerable to fluctuations in 
the projected galaxy density caused by foreground and/or 
background structures. 



Table [2] presents the results from letting 
100" < r min < 150" for Abell 1351 (with the KS93 
K map peak as cluster centre). It is seen that whilst c 
is decreasing with smaller r m i n , r2oo and A/200 remain 



stable for different 



Also worth noticing is that 



for 



< 130", c becomes constrained. However, as 



Train > r s for the different starting radii, we cannot 
obtain further conclusions from these results. As r s is 
even smaller for Abell 1995, we did not repeat this test 
for the cluster. iDietrich et al.l (|2005f l experience similar 
problems when attempting to determine the concentra- 
tion parameter for Abell 222 and Abell 223, concluding 
that obtaining a reliable c from weak lensing data only is 
difficult, if not impossible. 



14 



K. Holhjem et al.: Weak lensing analysis of Abell 1351 and Abell 1995 



Table 2: Results from varying the inner radius from where the shear values of Abell 1351 are measured. 



f min 


r s 


C 


^200 


M 200 


No. of 


X VDof 


["] 


["] 




lAroMpc] 


[lO M fc^M ] 


galaxies 




100 
110 


76 ±* 

61 +48 
u -29 

54 +60 

°^ -26 
gg+83 
uu -26 

07 +°0 

°' -20 
oo +oo 
J0 -15 


4 q+3 6 

-2.2 

6 0+ 4 ' 7 

" u -2.8 


1 7R+U.U57 
i ' ,u -0.13 

1.71 ± 0.10 


8.5 +*« 
7.9 ± 1.4 


15 630 
15 582 


1.11 
1.45 


120 
130 
140 
150 


6-8 + 3 3 

6 2 +s:s 

°- Z -2.6 
H-2 


i 71 +0.10 
' -0.086 
1 7f; +0.057 

± -' u -0.11 
i 170+0.11 
1 -'°-0.13 

1 74+o-ii 
-0.10 


79+1.4 
' - y -1.2 

8.5+™ 

"is 

8-311-2 


15 529 
15 482 
15 428 
15 358 


1.31 
1.33 
1.16 
0.96 



6.2.1. Best fit concentration parameter 



Bullock et alJ l|2001l ) present dark matter halo simula- 
tions, attempting to find a "best fit concentration param- 
eter" applicable to all types of halos. They find that for 
halos of the same mass, the concentration, c V i r = r vil /r s , 
can be given by c V h- oc (1 + Zd) _1 - This is contrary to ear- 
lier beliefs that c v - lr does not vary much with redshift. 
Numerically simulated massive clusters typically hav e con- 
centration parameters ^4 — 5 ([Bullock et al.ll2001r ). This 
is within the limiting values for both clusters, although 
looking at Fig. [TJ the outcome from varying c seems to 
better follow the shear values of Abell 1351. 

There exists several examples of large concentra- 



tion 
find 



parameters in the literature. iKneib et al 



22 



+9 



for the central mass 



concentration of 



the cluster CI 0024 + 1654. iGavazzil (1200511 conclude on 



11.73 ± 0.55 for MS2137 - 23, while 



Broadhurst et al 



2005) find c=13.7lj-? for Abell 1689. 



Limousin et al 



20061 1 present a thorough discussion of the different con- 



centration parameters derived for Abell 1689 in the litera- 
ture, and conclude that a distribution of best fit c param- 
eters is needed for observed lensing clusters in order to 
provide a sample large enough to make an adequate com- 
parison with simulations. A rece nt study of observed con- 
centra tion values for clusters by lComerford fe Natarajan 
1 20071 1 show that the best fit lensing-derived c parameters 
are systematically higher than concentrations derived via 
X-ray measurements, a difference which can be at least 
partly explained by effects of triaxia li ty of cluster halos 
ilCorless fc Kind |2007l : IGavazzil 120051 : lOguri et alJ 120051 
Clo we et al.l 2004 ) or t he sub-structure within the clusters 
l|King fc Corlessll2007t h although the latter effect may also 
produce a negative bias of c values. In addition, baryonic 
physics can increase the concentration parameter mildly 
by up to 10% as compared to dissipatio nless dark matte r 
in pure dark matter simulations (see e.g. Lin et al. 20061 ). 



6.3. Centre position 

In addition to the three centre positions tested in Sec t.El 
we com puted re maps with the method of Kaiser fc Squires] 
1 1993L KS93) and utilised the peak of this surface mass 
distribution as a fourth cluster centre. The KS93 method 
assumes that 7 = g, which is not a good approximation 
near the centres of massive systems. Therefore, in compar- 
ison with the other methods, it provides us with a refer- 



ence point as for how large a variation one can reasonably 
expect for the various centroiding methods. 

All centre positions obtained with the four methods lie 
within 1' and hence represent the errors expected when 
using the peak of a re map as cluster centre. As can be 
seen from TableQ] varying the centre position only slightly 
(maximum distance between cluster centres is < 55") can 
lead to different mass estimates. Although within error 
bars, the results from fitting NFW using a fixed c varies 
most. The NFW fitting of two parameters is more stable 
with a smaller spread in A/200 • This is also reflected in 
X 2 /Dof, as a value closer to 1 is a better fit. 

Worth noticing is the generally smaller differences be- 
tween the results of Abell 1995 as compared to those of 
Abell 1351. The concentration parameter seems also bet- 
ter constrained for Abell 1995, where we could not obtain 
an upper limit for c only in the case where the BCG was 
used as the centre reference. For Abell 1351 on the other 
hand an upper limit for c could not be obtained for any of 
the cluster centres chosen. This is consistent with the fact 
that Abell 1351 is not in dynamical equilibrium, lacking a 
well-defined cluster centre. The results obtained from fit- 
ting spherically symmetric models will hence depend on 
the centre chosen. 



6.3.1. Strong lensing features 

This picture of the centre of Abell 1351 is confirmed also 
by strong lensing effects, which is in general susceptive to 
sub-structures in clusters. For both clusters recent archival 
WFPC2 HST data exists, taken for an ongoing snapshot 
survey of X-ray luminous cluster (HST PID 11103, PI: H. 
Ebeling). The images are taken through the F606W filter 
totalling 1200s exposure time each. 

Considering both the colours and morphologies of 
galaxies in our V — I data and the morphologies in the 
HST images, there are at least half a dozen plausible arcs 
and arclets visible in each of the two clusters. The lens- 
ing pattern for Abell 1351 appears to be very complex 
and does not indicate a single, well-defined centre. This 
is supported by the presence of several elliptical galax- 
ies which are of similar brightness as the BCG. On the 
contrary, for Abell 1995 several arc(let)s are well aligned 
around the BCG, apart from three which are obviously 
associated with individual cluster galaxies. Thus adopting 
the BCG for Abell 1995 as the centre is justified. Strong 
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lensing will therefore not offer more constraints on the 
determination of the centre of mass than we already have. 

6.4. The mass reconstructions 

In Sect. [4] we presented the weak lensing reconstruction of 
the clusters' surface mass density, detecting the clusters 
on the 5cr level. In case of Abell 1351 a significant neigh- 
bouring peak A was detected, which coincides with the 
galaxy distribution. 

We note sig nificant differences comparing our mass 
maps to those of lDahle etaH d2002h. who used the KS93 
inversion method ([Kaiser fe Squiredll993l ). Abell 1995 ap- 
pears rather circular in both reconstructions, with at- 
tached filamentary structures of low significance seen in 
the KS93 map. In their re-anal y sis of the Abell 901 su- 
percluster field iHevmans et al. 1 2008^ argue that such 
filamentary structures could be enhanced by the KS93 
method itself. This algorithm assumes 3 = 7 near the 
critical cluster cores, which hence overestimates n, and the 
smoothing implemented could then lead to the apparent 
merging of neighbouring peaks. 

Our map of Abell 1351 appears roughly circular at the 
position of the main cluster with a significant extension 
towards the north-east, wherea s it shows an exten sion to 
the south-west in the work of iDahle et ail <|2002h . Both 
reconstructions have peak B in common. Since it is not 
associated with any over-density of galaxies but appears in 
both reconstructions based on very different data sets, the 
most likely explanation is a chance alignment of galaxies 
triggering this detection. 

7. Conclusions 

Utilising observations from CFH12K we find the masses of 
Abell 1351 (Abell 1995) to be Af 200 = M(r 20 o) ~ 8 - 9 x 
10 14 h^M® (M 200 - 5 - 6 x 10 14 h^M Q ). These results 
are derived from fitting our data to the NFW profile, alter- 
ing both c and r 20 o, to find the best theoretical fit to our 
shear measurements (Fig. [7]). Our k reconstructions yield 
mass estimates of M 2D (< ^200) = H.7±3.1 x 10 14 h^M Q 
(M 2D (< 7-200) = 10.5 ± 2.7 x lO u h^M ) for Abell 1351 
(Abell 1995). 

Our results illustrate that using solely weak lensing 
measurements, with no photometric or spectroscopic in- 
formation to separate cluster members from background 
galaxies, the concentration parameter found for a galaxy 
cluster is poorly constrained. Future spectroscopic redshift 
measurements and strong lens modelling of the arcs seen 
towards the cores of Abell 1351 and Abell 1995 might 
help improve the constraints on their concentration pa- 
rameters. However, the change in c value resulting from 
small variations of the centre position of Abell 1351 indi- 
cates that the mass distribution in the core of this dynam- 
ically unrelaxed cluster is too complex to be meaningfully 
fit by the NFW mass density profile. M 20 o, on the other 
hand, varies only slightly for both clusters when either 
r m ; n is decreased or the cluster centre varied. Although 



the KSB+ PSF correction leads to an underestimation of 
the cluster masses (with a maximum of 20% in our case), 
increasing the background galaxy ellipticities correspond- 
ingly still yields masses within error bars of our current 
estimates. We therefore conclude that the mass estimates 
presented are robust. 
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